function summation_of_amplitude_odn(CONSTS, plot_data)
    
    phi_vec = (-pi : pi/1000 : pi)';

    I_odn_pl_vec = sum_I_odn_pl(phi_vec, CONSTS);
    
    ang = angle(I_odn_pl_vec); % 1i*log(I_odn_pl_vec./(abs(I_odn_pl_vec)));
    ang_correct = ang;
    ang_correct(ang < -0.5*pi) = ang(ang < -0.5*pi) + 2*pi;

    if(plot_data)
        figure; plot(phi_vec*180/pi, abs(I_odn_pl_vec), 'b-'); 
        title('Abs of I_{\Sigma}(\phi) for odn plazma'); grid;
        xlabel('\phi, degrees'); ylabel('|I_{\Sigma}(\phi)|');

        figure; plot(phi_vec*180/pi, ang_correct*180/pi, 'b-'); 
        title('Angle of I_{\Sigma}(\phi) for odn plazma'); grid;
        xlabel('\phi, degrees'); ylabel('arg(I_{\Sigma}(\phi))'); 
    end

end